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We introduce a purely graph-theoretical object, namely the coding clique, to construct quantum error- 
correcting codes. Almost all quantum codes constructed so far are stabilizer (additive) codes and the construc- 
tion of nonadditive codes, which are potentially more efficient, is not as well understood as that of stabilizer 
codes. Our graphical approach provides a unified and classical way to construct both stabilizer and nonadditive 
codes. In particular we have explicitly constructed the optimal ((10,24,3)) code and a family of 1-error detect- 
ing nonadditive codes with the highest encoding rate so far. In the case of stabilizer codes a thorough search 
becomes tangible and we have classified all the extremal stabilizer codes up to 8 qubits. 



Introduction 

Quantum error correcting code (QECC) [2] [3] has be- 
come an indispensable element in many quantum informa- 
tional tasks such as the fault-tolerant quantum computation 
[5 1, the quantum key distributions [6], and the entanglement 
purification [7] and so on, to fight the noises. Roughly speak- 
ing, a QECC is a subspace of the Hilbert space of a system 
of many qubits (physical qubits) with the property that the 
quantum data (logical qubits) encoded in this subspace can be 
recovered faithfully, even though some physical qubits may 
suffer arbitrary errors, by a measurement followed by suitable 
unitary transformations. 

Almost all QECCs constructed so far are stabilizer or ad- 
ditive codes |9] [10 1 whose code space is specified by the 
joint +1 eigenspace of a stabilizer, an Abelian group of tensor 
products of Pauli operators. Special examples include the CSS 
codes [11 1, the topological codes 11121 . color codes IfTJI . and 
the recently introduced entanglement-assisted codes [14|. Be- 
cause of its structure a stabilizer code has always a dimension 
that is a power of 2 and is usually denoted as [[n, k,d]], which 
is a 2 k dimensional subspace of an n-qubit Hilbert space. Here 
d denotes the distance of a quantum code meaning that [^=i]- 
qubit errors can be corrected and all the errors that acts non- 
trivially on less than d qubits either can be detected or sta- 
bilize the code subspace. Apart from a classification of the 
[[«, 0, d]] codes up to 12 qubits [15|, the classification of sta- 
bilizer codes encoding a nonzero number of logical qubits has 
not been achieved yet. 

The code without a stabilizer structure is called as a non- 
additive code. Usually we denote by ((«, K, d)) a nonadditive 
QECC of distance d that is a ^-dimensional subspace of an 
n-qubit Hilbert space correcting [^]-qubit. Since less struc- 
tured than the additive codes, the nonadditive codes may be 
more efficient in the sense of a larger code subspace on one 
hand, harder to be constructed on the other hand. The first 
nonadditive code [16| that outperforms the stabilizer codes is 



the 1-error detecting code ((5,6,2)). Its generalization 1 17 1 is 
outperformed by another family of codes of distance 2 ifTHl . 
There are also some efforts to construct nonadditive codes 
systematically |[T9l l20l |2TI . Only recently the first error- 
correcting code, namely a ((9,12,3)) code, that outperforms 
the optimal stabilizer code is explicitly constructed in 11221 . 
Later on some other nonadditive codes have also been found 
HH. 

In this article we present a graphical approach to construct 
both additive and nonadditive codes in a unified and system- 
atic manner. Graph is related intimately to the construction 
of classical error-correcting codes (e.g. Tanner graph |24|) 
and also finds applications in the quantum stabilizer code via 
graph states E5]|26]|27l|28). with the hel P of g ra P h states, we 
relate a pure graph-theoretic object — the coding cliques — 
to the construction of quantum codes, especially of nonaddi- 
tive codes. An algorithm to search systematically for quantum 
codes via coding cliques is outlined. 



Graph and graph-state basis 

A graph G = ( V, T) is composed of a set V of n vertices and 
a set of edges specified by the adjacency matrix T, which is 
an nxn symmetric matrix with vanishing diagonal entries and 
T a b — 1 if vertices a, b are connected and T a b = otherwise. 
Here we consider only undirected simple graphs. The neigh- 
borhood of a vertex a is denoted by N a = {v e V\Y av = 1), i.e, 
the set of all the vertices that are connected to a. 

The graph states l25l l26l are useful multipartite entangled 
states that are essential resources for the one-way computing 
E9l and can be experimentally demonstrated l30l . Also the 
graph state plays the key role in our graphical construction of 
QECC. Consider a system of n qubits that are labeled by those 
n vertices in V and denote by X a , i/ a ,Za, and I a three Pauli 
operators and identity matrix acting on the qubit a e V. The 
«-qubit graph state associated with G reads 12511261 
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where \ft) z is the joint eigenstate of Z. a (a e V) with {—Yf a 
as eigenvalues, \+) v x is the joint +1 eigenstate of X a {a e V), 
and Uai, is the controlled phase gate between qubits a and b. 
The graph-state basis of the n-qubit Hilbert space refers 
to {|r c > := ZdO) with C e 2 V , the set of all the subsets 
of V. A collection of subsets \C\, C2, ■ ■ ■ , Ck) specifies a K 
dimensional subspace of ti n that is spanned by the graph-state 
basis {ITc,)}^! ■ The graph state |T) is also the unique joint +1 
eigenstate of n vertex stabilizers 

Ga=XaY[Zb-=XaZN a , a € V. (2) 

beN a 

Obviously a \r c ) = -|T C > if a e C and g a \T c ) = \T C ) other- 
wise. For a given subset S Q V the operator Q$ '■- Ylaes Qa 
also stabilizes the graph state, i.e., QsW) = |F). 

A graph determines uniquely a graph state and two graph 
states determined by two graphs are equivalent up to some lo- 
cal Clifford transformations (LCTs) iff these two graphs are 
related to each other by local complements (LCs) I5TI . An 
LCT is a special local unitary transformation that maps Pauli 
operators to themselves up to some sign changes. An LC of 
a graph on a vertex v refers to the operation that in the neigh- 
borhood of v we connect all the disconnected vertices and dis- 
connect all the connected vertices. All the graphs on up to 
12 vertices have been classified under LCs and graph isomor- 
phisms IT31 . 

Coding cliques and QECC 

From the standard theory of the QECC we know that if a 
set of Pauli errors can be corrected then all the errors that are 
linear combinations of the Pauli errors belonging to that set 
can also be corrected. One crucial observation is that when 
acting on the graph states, because of the vertex stabilizers in 
Eq.(|2]), the Pauli errors can be equivalently, up to some signs, 
replaced by phase flip errors, which motivates us to introduce 
the notion of coding cliques as follows. 

Given a graph G = (V, T) we define the neighborhood of 
a subset S c V to be Ns = A ve sN v , where we have denoted 
by C A D := C U D - C n D the symmetric difference of two 
sets C,D c V. Given any integer d (1 < d < ri) we define a 
d-purity set as 

$4 = Is QV\\S UN s \<d\. (3) 

Here \S\ denotes the number of the elements in set S . For a 
subset of vertices C c V if there exist two subsets a>,6 c V 
such that C = 6 A N w and \a> U 5\ - p, we call the ordered pair 
(6, cj) as a p-cover of C or equivalently C is said to be covered 
by p vertices or by (6, to). We denote by B> d the d-uncoverable 
set that contains all the subsets of V which cannot be covered 
by less than d vertices, i.e., 

B d = 2 v -ldANj\ \(0U6\ <d). (4) 



A coding clique of a graph G is a collection of vertex sub- 
sets \C\ , C2, . . . , Ck] that satisfies 

Condition 0. € Cj; 

Condition 1. |C; n S\ is even for all 1 < i < K and S e S d ; 

Condition 2. C, A Cj e B d for all 1 < 1 \ + j < K. 

A ^-clique of a given graph refers to a subset of K vertices 
that are pairwise connected. Our coding clique C? of a graph 
G is exactly a /T-clique of the super graph G defined as fol- 
lows. The vertices of the super graph G include the empty 
set and all the nonempty vertex subsets that belong to and 
satisfy Condition 1. Two different vertices C,C of the su- 
per graph G are connected if C A C e D^. Thus the coding 
clique is a pure graph-theoretic object that is formulated 
in a constructive way. The following theorem relates the cod- 
ing cliques to the construction of QECCs. (See Appendix for 
proofs.) 

Theorem 1 Given a coding clique of a graph G on 
n vertices, the subspace spanned by the graph-state basis 
{\T C )\C 6 Cf}, denoted as (G, K, d), is an ((n, K, d)) code. 

As an application, a systematic search for the quantum 
codes can be done according to the following algorithm: i) 
To input a graph G = (V, F) on n vertices that may be con- 
nected or disconnected; ii) To choose a distance d and com- 
pute the of-purity set S d and the of-uncoverable set so that 
a super graph G can be built; iii) To find all the ^-clique 
of the super graph G; iv) For every coding clique we obtain a 
(G, K, d) code , i.e., an ((n, K, d)) code that is spanned by the 
graph-state basis {|r c > = ZcW)\C e Cf}. 

It is not difficult to see that if the c/-purity set S d is empty 
then the (G, K, d) code is pure, i.e., every error acting nontriv- 
ially on less than d qubits can be detected. 

If a coding clique form a group with respect to the sym- 
metric difference A, then this coding clique is referred to as 
a coding group of the graph. Because the coding group is 
an Abelian group with self inverse, the number K of its ele- 
ments must be a power of 2, i.e., K — 2 k for some integer k, 
which is referred to as the dimension of the coding group. A 
^-dimensional coding group has k independent generators. If 
we find a coding group then we obtain a stabilizer code and 
all the stabilizer codes can be found this way because of the 
following theorem. 

Theorem 2 Every coding group C 2 d (G) provides a stabi- 
lizer code [[n, k, d]], denoted as a [G, k, d] code. Every stabi- 
lizer code [[n, k, d]] is equivalent to a [G, k, d] code for some 
graph G on n vertices. 

From the same graph we may obtain inequivalent (G, K, d) 
codes and different graphs may provide equivalent codes. To 
reduce the number of the graphs to be searched we have also 
investigated how the coding clique changes under the local 
complements of the graph. For convenience we denote by G v 
the graph obtained from G by making an LC on vertex v. For 
a given vertex v and a subset C c V we denote C v = C if v ^ C 
and C v = C A N„ if v e C. 
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LC rule for coding cliques If Cf is a coding clique (group) 

of the graph G, then Cj = {C V ,|C e Cj] is the coding clique 
(group) of the graph G v . Two codes specified by coding 
cliques (groups) and ofG and G v are equivalent under 
LCTs. 

Two quantum codes are regarded to be equivalent if they 
are related to each other by LCTs plus permutations of physi- 
cal qubits. Therefore we need only to take those inequivalent 
classes of graphs under LCs and graph isomorphisms as inputs 
to our algorithm. In what follows we shall document some re- 
sults obtained via the algorithm outlined above ll32l . A QECC 
will be specified by a graph together with a coding clique. For 
a stabilizer code we will only specify the generators of the 
coding group and for a nonadditive code we will list all the 
members of the coding clique. 

Nonadditive codes 

At first let us reproduce some known nonadditive codes via 
our graphical construction. The first example is the first non- 
additive 1 -error detecting code ((5, 6, 2)) [ 16 1 that outperforms 
the optimal stabilizer code [[5,2,2]]. The graph provides the 
code is the loop graph L 5 on 5 vertices as shown in Fig.l.(a) 
and the coding clique is {C,}^ =1 where 

d = 0, C 2 = (235), C 3 = {341), 
C 4 = (452), C 5 = (513), C 6 = (124). (5) 

The second example comes from the first 1-error correcting 
nonadditive code ((9, 12, 3)) that outperforms the optimal sta- 
bilizer code [[9, 3, 3]]. The code is specified by the loop graph 
Lg on 9 vertices and the coding clique contains those 12 sub- 
sets { V,]'f j defined as in Our search results show that the 
code (Ls,6,2) and the code (Lg, 12,3) are unique. In addition 
there is no (G, 13, 3) code for any graph G on 9 vertices. 

Another example is Rain's ((2m + 3, 6 • 4 m_1 , 2)) code 07), 
which can be specified by the graph made up of L5 and 2(m— 1 ) 
pairwise connected vertices as shown in Fig. 1(b). The coding 
clique contains all the subsets of form {C,- A t/]f =1 where C,'s 
are given in Eq.{5]) and U can be any one of 2 2m ~ 2 subsets of 
vertices generated by {2, aj} and {5, b[} for j, I = 1, 2, . . . , m— 1 
with respect to the symmetric difference A. 

Now let us construct some new codes. The first result is a 
family of 1-error detecting codes ((An + \,M„ + 1,2)) where 
M n = 2 4 " -1 - 5C 2 ". A family of nonadditive codes of dis- 
tance 2 [ 18 1 is constructed which encodes an M„ -dimensional 
subspace if there are An + 1 physical qubits. The code can be 
specified by the star graph on An + 1 vertices V"4„+i centered on 
vertex o as shown in Fig.l (c) (indicated by red edges). The 
coding clique reads 

Cc V 4 „ +1 |o i C;\C\ =2/, or|C|=2n + 2/+l,0</<n-l 

(6) 

For the star graph the dimension M„ of the code space is opti- 
mal 11181 . We consider instead the graph as shown in Fig. 1(c) 




(a) (c) (d) 



FIG. 1: The graphs and coding cliques for nonadditive codes, (a) Six 
elements of the coding clique for the ((5, 6, 2)) code, (b) A graph for 
Rain's ((2m + 3,6- A"'- [ , 2)) code, (c) The graph for the ((An + 1, M„ + 
1,2)) code together with an element of the coding clique, (d) The 
((10, 24, 3)) code. 



that is built on a star graph with red edges and some addi- 
tional black edges. The coding clique, in addition to Eq.([6]i, 
contains one more subset of all 2n + 1 blue-colored vertices 
{o,4Z - 2,4/ - 1}" =1 as indicated in Fig. 1(c). Therefore we 
obtain a family of nonadditive code of distance 2 with an en- 
coding rate M„ + 1 . 

The second result is the optimal ((10, 24, 3)) code. The lin- 
ear programming bound |10| indicates that the largest code 
subspace of a 1-error correcting code on 10 qubits is of di- 
mension 24. Among 3132 inequivalent connected graphs and 
all the disconnected graphs on 10 vertices there is a unique 
(G10, 24, 3) code as depicted in Fig. 1(d) in which blue-colored 
vertex sets indicate all the nonempty elements of the coding 
clique. 



Stabilizer codes 

Because of Theorem 2 the search on all the inequivalent 
graphs on n vertices for all the coding groups according to our 
algorithm provided above will exhaust all the stabilizer codes 
of n qubits. Therefore a classification of all stabilizer codes 
is tangible. We consider in the following the classification 
of those extremal codes, i.e., the code with a maximal k or a 
maximal d given n, up to 8 qubits, as listed in the table in iflOl . 

Up to six qubits, there are three extremal codes [[4,2,2]], 
[[5,1,3]], and [[6,1,3]], which are all unique and corre- 
sponding graphs and coding groups are shown in Fig. 2. 
For [[4,2,2]] we have shown 4 different graphs and coding 
groups, all resulting an equivalent code. There is also a trivial 
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[[5,1,3]] [[6,0,4]] [[8,1,3]] [[10,1,4]] [[10,2,4]] 




[[S,3,3]]„ [[S,3,3]] x [[8,3,3]],, [[8,3,3]] # [[8,3,3]]_ [[S,3,3]] / 



FIG. 2: The classification of stabilizer codes part I. There are three 
primary colored sets of vertices: red set R, yellow set Y, and blue set 
B. Orange, green, purple, and black vertices belong to R n Y, Y n B, 
BCiR, and RnYnB, respectively. Each primary color set represents 
a generator of the coding group. 



[[7, 1,3]] codes. 

Among 101 inequivalent connected graphs on 8 ver- 
tices there are only 6 LU-inequivalent graphs that admit 3- 
dimensional coding groups as shown in Fig. 2. In each graph 
there may be many different coding groups and only one cod- 
ing group is shown for each graph. Disconnected graphs 
on 8 vertices provides no further [[8,3,3]] codes. From the 
weight distributions of corresponding codes it can be proved 
that the code [[8, 3, 3]] is unique (Appendix). In Fig. 2 we also 
show three different graphs and the coding groups that pro- 
vide codes that are equivalent to Shor's [[9,1,3]] code (T) 
and a LC-equivalent graph for the impure [[10,1,4]] code 
constructed in |26l . In comparison we have found a pure 
[[10, 2,4]] code with a weight distribution (90 6 , 135 8 , 30 K) ). 



Discussion 




(2) * (4) (25) 






OOO 
(351 (37) 





FIG. 3: Classification of stabilizer codes part II. The table of 16 in- 
equivalent [[7, 1,3]] codes with 10 different weight distributions sep- 
arated by lines and dashed lines. Graphs with blue vertices provide 
pure codes and graphs with red vertices provide impure codes. 



[[6, 1, 3]] code obtained from [[5, 1, 3]] code by appending an 
isolated vertex to the pentagon as noticed in [ 10 1 . In the clas- 
sification of 7-qubit codes we will neglect also those codes 
obtained from [[5, 1,3]] and [[6, 1,3]] by appending a discon- 
nected subgraph. 

Among 26 inequivalent connected graphs and on 7 vertices 
there are 14 graphs that admit a 1 -dimensional coding group 
as shown in Fig. 3, where the numberings of the inequivalent 
graphs are adopted from [25|. One graph may provide in- 
equivalent codes and different graphs may provide equivalent 
codes. Steane's 7-qubit code [4] is provided by graph No.43 
with the weight distribution Wg. (A detailed account of weight 
distribution [33 1 is given in Appendix.) It is interesting to no- 
tice that a [[7, 1, 3]] code can be provided by a disconnected 
graph, whose two subgraphs cannot admit any coding clique 
alone. Disconnected graphs on 7 vertices provides no further 



We provide a classical recipe to cook quantum codes, which 
is exhaustive for stabilizer codes and systematic for nonaddi- 
tive codes. A classification of stabilizer codes has been done 
up to 8 qubits. Since all the known good nonadditive codes 
can be reproduced and an optimal code ((10,24,3)) is con- 
structed via our graphical approach, we speculate that our al- 
gorithm is also exhaustive for nonadditive codes. 

With an improved computational power a classification of 
the stabilizer codes on more qubits and the discovery of more 
good codes either stabilizer or nonadditive codes are expected. 
It should admitted that the clique finding problem is intrinsi- 
cally a NP-complete problem. Therefore analytical construc- 
tions of the coding cliques, at least for some special family 
of graphs such as loop graphs and hyper cubic graphs, de- 
serve exploring. A direct construction of a series nonadditive 
1 -error detecting codes with highest encoding rate so far de- 
scribed here may provide a clue. In addition it is not difficult 
to generalize the idea of coding cliques to the construction of 
nonbinary codes via graph states for systems with more than 
2 levels. 

Our precious quantum data can be protected from decoher- 
ences either in a dynamical manner as in the decohere nce-free 
subspace approach, or in a geometric manner as in geometric 
computations, or in a topological manner as in the topological 
quantum computations based on the topological codes, a spe- 
cial stabilizer codes. Our results make a bridge between the 
exciting classical field of graph theory and the quantum error 
correction. Equipped with the one-way computation model 
based on graph states and the graphical QECCs, we may en- 
vision a graphical quantum computation based directly on the 
graphical objects. 

SXY acknowledges financial support from NNSF of China 
(Grant No. 90303023 and Grant No. 10675107), CAS, and 
WBS (Project Account No): R-144-000-189-305, Quantum 
information and Storage (QIS). 
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Appendix 

Proof of Theorem 1 It is enough to prove that for any 
error acting nontrivially on less than d qubits we have 
<r c |£ d |r C '> = f(& d )5 ca miM for all C,C 6 C* With- 
out lose of generosity we assume that S d = X^Zs for some 
pair of subsets 6, to with \co U S\ < d, which represents that 
there are X, J/, and Z errors on the qubits in to - 5 Pi to, to Pi 5, 
and 6 - to n 5 respectively. When acting on the graph state 
the error & d oc Q^Zsi can be replaced by phase flip errors Zsi 
on Q := 5 A N w . If Q is empty then 6 = N w and the error is 
proportional to Q w . In this case we have QJ\Tc) = Wc) for 
all C G C k d because \a> n C\ is even which stems from the fact 



that \to U NJ < d, i.e., oj G 



and Condition 1. Thus the 



error behaves like a constant operator on the coding subspace 
and can be neglected. If Q. is not empty then O ^ Dj U 
because it is covered by (5, to) and \6 U a>\ < d. As a result 
<r c |£rf|r c <> oc (riZcAcZnlO = for all C,C e C* U be- 
cause condition 2 ensures that C AC + O. 

Proof of Theorem 2 We prove the second part first, i.e. to 
construct a graph and a coding clique from a given stabilizer 
code. We do not need to introduce those input vertices as in 
Ref.GS). 

A stabilizer is a set of commuting observables that are 
tensor products of Pauli operators \X, J/, Z) and the iden- 
tity operator I on each qubit. A stabilizer code [[«, k, d]] 
is the simultaneous +1 eigenspace of n — k generators 
{S\,S2, ■ ■ ■ ,S n -k) of the stabilizer. Stabilizer codes that can 
be related to each other by a local Clifford transformation 
(LCTs) and permutations of physical qubits are equivalent. 

With the help of the binary representation (I — > 00, X — > 
10, Z — > 01, J/ — > 11) of the local operators we can describe 
the stabilizer a check matrix [G X |G ; ] where (n - k) x n ma- 
trix G X ( Z ) is formed by all the first (second) digits of the local 
operators of the generators of the stabilizer. By relabling or a 
different choice of the generators of the stabilizer, permutating 
the qubits, and making suitable LCTs, one can always bring 
the check matrix into its standard form f34ll 



[I r ,A rxk \D rxr + AE' , E rxk ] (r = n-k) 



(7) 



where denotes the identity matrix and the submatrix D is 
symmetric as required by the commutativity among the gen- 
erators, i.e., G X G\ being symmetric. Here all the additions are 
additions module 2. 

A basis of the code subspace can be constructed by adding 
to the stabilizer additional k independent generators. The most 
general choice, up to an LCT and a choice of different gener- 
ators, is to add the k generators specified by the check matrix 

[0J k \E T + FA T ,F kxk ] (8) 

where F is an arbitrary symmetric matrix. By adding the 
check matrix Eq.(|8]l multiplied by A from left to the the stan- 
dard form Eq.([7]i we obtain the check matrix of n generators 

as 



D + AFA T 
E T + FA T 



E+AF 
F 



If one of the diagonal element F, say T mm , is nonzero then 
we can make it vanish by performing a LCT (if — > X — > 
-J/; -Z — > -3) to the m-th qubit. Therefore we obtain a stan- 
dard adjacency matrix F of a graph: n x n symmetric matrix 
with vanishing diagonal entries. In terms of the vertex stabi- 
lizers \Q a } of this graph, the stabilizer is generated by 



Qa Y^Qn-k+ct" a= l,. ..,n-k. 



(10) 



Denote by C„, the set of vertices on which the entries of the 
m-th row (m = 1, ...,£) of the k x n matrix [A T , I k ] do not 
vanish and by C the set of subset of V generated by {C m } k m=l 
with respect to the symmetric difference A. Then the graph- 
state basis {Zc|r)|C G C) spans the [[«, k, d]] code. By def- 
inition C is a group with respect to the symmetric difference 
and G C. Secondly for every S e S,/ the graph stabilizer 
Qs , which can be a possible error, must stabilize the code sub- 
space, i.e., \S nC| must be even for all C e C. Thirdly, had one 
nonempty CeCa /?-cover (to, S) with p < d, the error opera- 
tor XojZ,s would be incorrigible because XaZslG) oc Zc\G). 
Thus C c Drf, i.e., Condition 3 is satisfied. As a result C is a 
coding group. 

One the other hand if we have a k dimensional coding group 
of a graph G which is generated by (C\, C2, ■ ■ ■ , C k ), then we 
have a code (G,2 k ,d) according to Theorem 1. Since k con- 
straints \S U Ci\ =even f or i = 1,2, ... ,k have exactly n - k in- 
dependent solutions which are denoted as (Si, S2, • • • <S n -k), 
we obtain a stabilizer of the code that is generated by {QstYlZl- 

Proof of LC rule for coding cliques Given a graph G — 
(V,T) and a coding clique we denote = {C V \C e Cj) 
with C v — C if v t C and C v — C A N v otherwise where N v 
is the neighborhood of v in the graph G. For convenience we 
denote A A'" B = A if m is even and A A m B = A A B if m is 
odd. Then we have concisely C v = C A' i nc ' N v , where we have 
denoted the single vertex set containing v also by v. After 
a local complement on vertex v made to G we obtain a new 
graph G v = (V, T). The neighborhood of a vertex set S c V in 
the new graph G v reads 



N s = N s A (5 n N v ) A lsnN ' ] N v 



(11) 



:= ii n \n 



(9) 



where Ns is the neighborhood of S in the original graph G. 

Firstly since G Cj and v ( we have e Cf . Condition 
is satisfied. Secondly, by denoting S v = S A' v ' v we see 
that if S G Sj(G v ), i.e., S belongs to the cZ-purity set of G v then 
S belongs to the cZ-purity set of G, i.e., S v 6 Ed(G), because 
S V U N Sr = S U N s ■ As a result \S „ n C| must be even for all 
C G Cf and S £ S d (G v ). Because 

(S n C v ) a (S v n C) = A |5nAg (v n C) A |vnc| (S n N v ) (12) 

is an even set, it follows that |S n C v | must also be even for all 
C v G C%. Thus Condition 1 is satisfied. Thirdly let us suppose 
that C v A C' v has a p-cover (6, to) in the new graph G v with 
p < d, i.e. C v A C' v = 8 A N w . Then (£', a/) provides a cover 
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for C AC" in G, i.e. C AC' = 6' A where 6' = 5 A (w n iV v ) 
and w' = w A™ v where 

m = \conN v \ + \vn(CAC')\ (13) 

Sincem + |vn5| is even we obtain |(5'Ua/| = p. ThusCAC'has 
a p < d cover in graph G which is in contradiction with Condi- 
tion 2 of a coding clique. Therefore C v A C' v cannot be covered 
by less than d vertices and belongs to the af-uncoverable set of 
G„, i.e., C v A C' v e B d (G v ) for every C,„ CJ 6 Cf (G). Condi- 
tion 2 is also satisfied. In addition if is a coding group of 
G then Cj is a coding group of G v since C v A C' v = (C A C') v . 

Let us denote by P and P two projectors of the code sub- 
spaces specified by C K d of G = (V,F) and of G v = (V,f) 
i.e., 

p= J] ZcirxriZc, p= J] Z c jr>(f|Z Cv . (14) 

CeC* C,,eCJ 

By denoting an LCT as U — y/-iX v YlueN,. yiZu we have 
'Z/|r><r|'Z/ t = |r>(F| and UZc^ = Zc = Zc, if v g C and 
WZcU 1 = iZc r @v if v e C. As a result Wl/ 1 = P, i.e., two 
codes are related to each other by the LCT 1A. 

Weight distribution and Classification of all the [[7, 1,3]] 
codes Given a K dimensional subspace with projector de- 
noted by P of n qubit Hilbert space and an arbitrary set u> 
of qubits, one can build an invariant [35| under local unitary 
transformations as 

Au = F Z |Tr(£P)|2 (15) 

supp(£)=w 

where the summation is taken over all Hermitian Pauli er- 
rors that acting nontrivially on the qubits in to. If a> = 
we have Ao = 1. If we sum over all possible subsets con- 
taining d qubits then we obtain the weight distribution 11331 
(Ao, A\, Az, . . . , A„) of a code where 

A d =Yj A <^ d = 0,l,...,n (16) 

are invariant under LUTs and permutations of qubits. Ob- 
viously Y!d = o^-d = 2"/K and we neglect Ao and those zero 
entries sometimes in the weight distributions. For nonaddi- 
tive codes A^ may be a fractional. For example the weight 
distribution of the nonadditive (Gio, 24, 3) code is ((y )6, 35s) 
meaning that Ao = 1, Ag = 20/3, Ag = 35, and A,- = other- 
wise. 

For stabilizer codes Aj is an integer which equals to the 
number of the stabilizers with the same weight d. The weight 
of a Pauli operator is the number of qubits on which it acts 
nontrivially. For all the stabilizer codes [[7, 1,3]] obtained by 
searching the coding groups can have only 10 different weight 
distributions as documented in the following table 
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By using the weight distributions we can only identify 10 dif- 
ferent classes of [[7, 1, 3]] codes. Further classification is done 
by the frequency analysis described as follows. Given d and a 
subset S of qubits we define 

F d (S)= J] A - d=\,2,...,n (17) 

to be the frequency of S , which is obviously an LU-invariant 
quantity. For every d we order all the frequencies Fd(S) of 
S containing the same number of qubits by their magnitudes. 
The resulting ordered series of frequencies is invariant under 
permutations of qubits. Any difference between the corre- 
sponding series of two codes will witness their inequivalency. 
As a result of frequency analysis we obtain 16 different in- 
equivalent codes and within each equivalent class all codes 
can be related to each other via explicit local Clifford trans- 
formations. 

Proof of the uniqueness of the code [[8, 3, 3]] All the sta- 
bilizer codes specified by the coding groups of graphs on 8 
vertices have the same weight distribution (286, 3s) and 3 sta- 
bilizers with full support form a group. Therefore the code 
must be pure (Ai,A2 = 0) and it is possible to bring two 
weight 8 generators of the stabilizer to Xy and Zv by LCTs. 
{V denotes 8 qubits here.) As a result one can always bring the 
local operator acting on the first qubit of the remaining three 
generators to the identity operator by choosing different gen- 
erators. Denote [G V |G Z ] as the check matrix of the remaining 
3 generators. The 3x8 matrix G x or G z must have distinct 
columns because all single-qubit Z or X errors can be cor- 
rected and two matrices G x and G z cannot have 2 or more 
identical columns on the same qubit because all the single- 
qubit if errors can be corrected and the code is pure. Thus the 
columns of G, form a map of the the columns of G x with only 
one fix point and this map is unique up to LCTs and permuta- 
tions | fOl . As a result the code [[8, 3, 3]] is unique. 



[1] P.W. Shor, Phys. Rev. A 52, R2493 (1995). 

[2] C.H. Bennett, D.P. DiVincenco, J.A. Smolin, and W.K. Woot- 

ters, Phys. Rev. A 54, 3824 (1996). 
[3] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 
[4] A. Steane, Phys. Rev. Lett. 77, 793 (1996). 
[5] E. Knill, R. Laflamme, W. H. Zurek, Science 279, 342 (1998); 

D. Gottesman, Phys. Rev. A 57, 127 (1998). 
[6] C. H. Bennett, G.Brassard, Proceedings of IEEE International 

Conference on Computers, Systems, and Signal Processing, 175 

(1984); A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 



7 



[7] S. Glancy, E. Knill, and H. M. Vasconcelos, Phys. Rev. A 74, 

032319 (2006). 
[8] D. Gottesman, Phys. Rev. A 54, 1862 (1996). 
[9] A. R. Calderbank, E. M. Rains, P. W. Shor, and N. J. A. Sloane, 
Phys. Rev. Lett. 78, 405 (1997). 
[10] A.R. Calderbank, E.M. Rains, P.W. Shor, and NJA. Sloane, 

IEEE Trans. Inf. Theory 44, 1369 (1998). 
[11] A. R. Calderbank and P.W. Shor, Phys. Rev. A 54, 1098 (1996); 

A.M. Steane, Phys. Rev. A 54, 4741 (1996). 
[12] A.Yu. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 
[13] H. Bombin and M. A. Martin-Delgado, Phys. Rev. Lett. 97, 
180501 (2006). 

[14] T. Brun, I. Devetak, M.-H. Hsieh, Science 314, 436 (2006). 
[15] L.E. Danielsen and M.G. Parker, J. Combinatorial Theory A 

113, 1351 (2006). 
[16] E.M. Rains, R. H. Hardin, P.W. Shor, and N.J.A. Sloane, Phys. 

Rev. Lett 79, 953 (1997). 
[17] E.M. Rains, IEEE Trans. Inf. Theory 45, 266 (1999). 
[18] JA. Smolin, G. Smith and S. Wehner, arXiv: quant- 

ph/0701065. 

[19] V.P Roychowdhury and E Vatan, Quantum computing and 

quantum communications 325, Lecture Notes in Comput. Sci. 

1509, (Springer, Berlin, 1999); arXiv: quant-ph/97 10031. 
[20] V. Arvind and K.R. Parthasarathy, arXiv: quant-ph/0206174; V. 

Arvind, P.P. Kurur, and K.R. Parthasarathy, quant-ph/02 10097. 
[21] V. Aggarwal and R. Calderbank, arXivxs/0610159. 
[22] S. Yu, Q. Chen, C.H. Lai, and C.H. Oh, arXiv: quant- 

ph/0704.2122. 

[23] A. Cross, G. Smith, J.A. Smolin, and B. Zeng, arXiv: quant- 
ph/0708.1021. 

[24] R.M. Tanner, IEEE Trans. Inf. Theory IT-27, 533 (1981). 
[25] M. Hein, J. Eisert, and H. J. Briegel, Phys. Rev. A 69, 062311 
(2004). 

[26] D. Schlingemann and R. E Werner, Phys. Rev. A 65, 012308 
(2002). 

[27] D. Schlingemann, Quant. Inform. Comput. 2, 307 (2002) 
[28] M. Grassl, A. Klappenecker, and M. Rotteler, IEEE Int. Symp. 

Inform. Theory Proceedings, 45 (2002). 
[29] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 5188 

(2001). 

[30] P. Walther, et.al., Nature 434, 169 (2005); C.Y Lu, et.al., Nature 

Physics 3, 91 (2007). 
[31] M. Van den Nest, J. Dehaene, and B.D. Moor, Phys. Rev. A 69, 

022316(2004). 

[32] We have used the clique finding program Cliquer to search for 
the coding cliques. Ref.: S. Niskanen and P.R.J. Ostergard, Cli- 
quer User's Guide, Version 1.0, (Communications Laboratory, 
Helsinki University of Technology, Espoo, Finland, Tech. Rep. 
T48, 2003). 

[33] P. Shor and R. Laflamme, Phys. Rev. Lett. 78, 1600 (1997); 

E.M. Rains, IEEE Trans. Inf. Theory 44, 1388 (1998). 
[34] R. Cleve, Phys. Rev. A 55, 4054 (1997). 

[35] M. Van den Nest, J. Dehaene, and B.D. Moor, Phys. Rev. A 70, 
032323 (2004). 



